PERPETUITIES WITH THIN TAILS, REVISITED 

PAWEL HITCZENKO t AND JACEK WESOLOWSKI 



f^"*) Abstract. We consider the tail behavior of random variables R which are solutions of 

t— I the distributional equation R = Q + MR, where (Q, M) is independent of R and \M\ < 1. 

Goldie and Grubel showed that the tails of R are no heavier than exponential and that 
if Q is bounded and M resembles near 1 the uniform distribution, then the tails of R are 
,-0 Poissonian. In this paper we further investigate the connection between the tails of R and 

P^P the behavior of M near 1. We focus on the special case when Q is constant and M is 

non-negative. 

OO 

i — i 1. Introduction 

In this note we consider a random variable R given by the solution of the stochastic 
ri equation 

(1.1) R = Q + MR, 

where (Q, M) are independent of R on the right-hand side. Under suitable assumptions 
^ on (Q, M) one can think of R as a limit in distribution of the following iterative scheme 

,> (1.2) R n = Q n + M n R n _ u n>l 

where R is arbitrary and (Q n ,M n ), n > 1, are i.i.d. copies of (Q,M), and (Q n ,M n ) 
is independent of R n -i- Writing out the above recurrence and renumbering the random 



(N 



variables (Q n , M n ) we see that R may also be defined by 



oo j—1 

(1.3) R^J2^U M ^ 

^ 3=1 k=l 

provided that the series above converges in distribution. Sufficient conditions for the almost 
sure convergence are known and have been given by Kesten [13] who also considered a 
multidimensional case when M is a matrix and Q a vector. For a nice detailed discussion 
of a one dimensional case we refer to the paper by Vervaat |19j : we only mention briefly 
here that Elog + \Q\ < oo and Elog \ M\ < suffice for the almost sure convergence of the 



series in (1.3) 



In the form (1.3) R has been studied in insurance mathematics under the name perpe- 



tuity. Since schemes like (1.2) are ubiquitous in many areas of applied mathematics, the 
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properties of R have attracted a considerable interest. We refer to [SI El El El E321 EIS1 EES] 
and references therein for more information and sample of applications. For examples of 
more recent work on perpetuities and their applications see [Tj EJ [HI [H] . A few additional 
situations in which perpetuities arise will be mentioned below. 

The main focus of research is the tail behavior of R. Kesten [13] showed that if P(|M| > 
1) > then R is always heavy-tailed. More precisely, he showed that if there exists a k 
such that E|M| K log + \M\ < oo, E\Q\ K < oo, and E|M| K = 1 then for some constant C 

F(\R\ > x) ~ Cx~ K , as x oo. 

Here, and throughout the paper the symbol f(x) ~ g(x) means that the ratio goes to 1 
as x — > oo. His result was rediscovered, reproved, and extended by several authors (see 
[3 [HI HQ]). In the complementary case, P(|M| < 1) the picture is much less clear. The main 
work we are aware of is that of Goldie and Griibel [8] who showed that in that case, the 
tails are never heavier than exponential and that if M behaves near 1 as a uniform random 
variable then the tails have Poissonian decay. In their arguments Goldie and Griibel relied 



on inductive arguments applied to (1.2) 



The main purpose of this note is to use systematically their approach to obtain additional 
information on the links between the behavior of M near 1 and the tail behavior of R. 
Following Goldie and Griibel (and also customs in large deviation theory) we will be 
interested in the asymptotics of the logarithm of the tail probability, i.e. lnP(|i?| > x) 
as x — > oo. Since we are mainly interested in establishing the links between M and R, 
we will often make additional, but common, assumptions when necessary. For example, 
we generally assume that Q and M are independent or even that Q = q is non-random. 
The independence assumption is typically needed only for the lower bounds on the log 
of the tail probability, the upper bounds are usually obtainable without it. Once the 
independence of Q and M is assumed the restriction that Q is degenerate does not seem to 
be a major restriction, but makes some of the arguments more transparent. It is rather the 
assumption that Q is bounded, which seems to play the more important role. Similarly, we 
will assume that M and q are non-negative. How the non-negative case differs from the 
general is relatively well understood (see e.g. arguments in [HI Theorem 2.1, Theorem 3.1, 
Lemma 5.3]) to see how arguments for non-negative case can be extended to more general 
situations. 

We would like to mention an interesting connection of perpetuities with a subclass of 
infinitely divisible laws, namely, as was shown by Jurek [12] all self-decomposable random 
variables (we refer to [T2] for the definition) can be represented as perpetuities R given 



by (1.1) with < M < 1. As a matter of fact, much more is shown in [T2], namely, if 



R is self-decomposable then for every random variable < M < 1 there exists a random 



variable Q (typically not bounded) such that (1.1) holds with (Q, M) independent of R on 



the right-hand side. This curious result seems to be of little help as far as general theory of 
perpetuities goes. In fact, one can take M to be any constant M — m G (0, 1) and equally 
well represent a self-decomposable random variable as a series of weighted i.i.d. random 
variables, with weights forming a geometric progression. Nonetheless, we mention that 
building on an earlier work of Thorin [T51 ITT] , Bondesson [3] proved a general result which 
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implies, in particular, that all gamma, inverse gamma, Pareto, log-normal, and Weilbull 
distributions are self-decomposable. Some of these results were obtained earlier by other 
authors and we refer to Bondesson [3], Section 5] for credits and more examples. 

2. General outline 

To begin the discussion, assume that \M\ < 1. Trivially, if \Q\ < q and \M\ is concen- 
trated on a proper subinterval (0, 1 — 5), 5 > of (0, 1) then the perpetuity R is a random 
variable whose absolute value is bounded by q/5 and thus has a trivial tail in the sense 
that P(|-R| > x) = for x > q/5. On the other hand if M is not bounded away from 1 
then we have the following observation due to Goldie and Griibel: 

Proposition 1. For 5 G (0, 1) let p$ := P(l — 5 < M < 1). Then for every such 5 and for 
all y > we have 

(2.1) F(R> q -(l-(l-5Y))>pl 

In particular, if for c G (0, 1) and x > q we set 

cq ln(l - c) 

o = — and y — 



(2.2) F(R >x)>[p^) in(1 - cqM = exp ( ^pi — ^ 



x " ln(l — cq/x) ' 

then we get that 

' ln(l - c) 
ln(l — cq/x) 

Proof: This was observed by Goldie-Griibel: For a given 5 > we let 

r = r s = inf{n > 1 : M n < 1 - 5}. 
Then by non-negativity and ( |1.3[ ), on {r > n} we have 

n 

J2 >£) g(l = | (1 -(!-*)»). 

fc=l 

Therefore, for all n > 1, 

P(-R > ^(1 - (1 - <^) n )) > V(M k > 1 - 5, 1 < fc < n) = pj -1 . 



Hence, 



P(i?> ?(!-(! -5) y )) >p|, for ally >0 
o 



which proves (2.1); (2.2) follows by a simple calculation. □ 

It is clear from the above proposition that if ps is strictly positive for every 5 > then 
the perpetuity R has non-trivial tails. It is then the behavior of M near 1 that determines 
the nature of the tails of R. It appears that essentials of such a behavior are shared by a 
class of equivalent distributions in the following sense. 
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Let n and v be probability distributions on [0, 1]. For any 5 G (0, 1) we denote fi$ = 
/i((l — 5, 1]) and us = z/((l — 5, 1]) We say that the distributions [i and v are equivalent at 
1 if 

3 e > and < d < D < oo such that 

(2.3) We(0,e]: d<^<D. 

As we mentioned earlier, our goal here is to shed some additional light on the relation- 
ship between the behavior of the distribution of M in the left neighborhood of 1 and the 
tails of R. To accomplish that we will develop in a systematic way the approach of Goldie 



and Griibel. For the upper bound this approach relies on iteration of (1.2) to get a uni- 
form upper bound on the moment generating function of R n for all n > 1 and then use 
exponentiation and Markov inequality to translate this bound into bounds on the tails. 
We will develop this in the next section, but to give a flavor of this argument we provide 



the following illustration: Consider (1.1) and assume that Q, M, and R on the right-hand 
side of (1.1) are independent (that is of course stronger than the usual assumption that 
(Q, M) are independent of R). Also, assume that < M < 1 and that m := EM < 1. To 
get an upper bound on the moment generating function Ee zR of R, the principle of what 
Goldie-Griibel did is the following: for n > 1 we have 

Ee zR " = Ee* ( - Qn+MnRn - l) = Ee zQ Ee zMRn ~ 1 < Ee zQ {l + mEfe*- 1 - 1)} , 

where in the last step we use the fact that for s > 

(2.4) Ee sM < Ee sBin(1 ' m) = 1 + m(e s - 1). 

To set up an induction we seek a function A(z) such that 

(i) Ee^- 1 < A(z), and 

(ii) Ee zQ {1 + m(A(z) - 1)} < A(z). 

Solving (ii) gives 

(1 -m)Ee zQ 
B(z) := ^ '- — p=- < A(z), 

for z such that mEe z ® < 1. Now, B(z) is recognized as the moment generating function of 

SfcLi Qk where N = Geom(l — m) and is independent of the sequence Qk, k > 1. So if we 
start with any R for which (i) holds with B(z) in place of A(z) then the induction goes 
through and, under a reasonably weak assumptions on Q, we get an exponential upper 

bound on the tail of R. In particular if we take Q = 1 and M — Bin(l, m) then R has 
moment generating function bounded by that of a geometric random variable and hence 
sub-exponential tails as was already shown by Goldie and Griibel. 

We mention briefly that the sums described by B(z) are yet another example of perpe- 
tuities. Sums like these are of interest in renewal theory and risk assessment, for example. 
They have been studied before, for instance in [U [20], under the name geometric convo- 
lutions and geometric random sums, respectively. We refer the interested reader there for 
more information and further references. 
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As for the lower bound, the best that is available at this point is argument based on 
Proposition [TJ Interestingly, this proposition provides a surprisingly good lower bound. By 
this we mean the fact that if the upper bound obtained by the above method is constructed 
carefully so as to be relatively tight, then one can usually obtain a lower bound of a similar 
strength from Proposition [T} This will be seen in several situations below. It is thus 
important to understand how to construct a tight upper bound. Although we do not have 
a general result to that effect, in the last section we will provide an argument in a particular 
example that provides a heuristic which should work well in other cases. 

The rest of the paper is organized as follows, in the next section we will discuss an upper 



bound and in particular, we will state an inequality (see (3.6) below) that is crucial for the 
inductive argument. In subsequent sections we will illustrate this with several examples. 
Those include beta(a, /3) densities, and what (for the lack of a better name) we call the 
generalized beta(l,/3) densities. The reason for considering beta distributions is that one 
might reasonably hope that they provide a natural parametrization of a behavior of M 
near 1, which could be translated to the tail behavior of R. This, however, is not the case, 
since as we will show all beta distributions lead to the same, namely Poissonian, behavior. 
It turns out that a much more rapid than power-type variability of M at 1 is needed 
to observe a different tail behavior of R. We will then construct densities for which the 
logarithm of the tail probability will have power behavior — x r , for 1 < r < oo. In the last 
section we will discuss one more example mainly to illustrate a techinque of constructing 
M that would give a particular tail behavior of R in other situations. 

3. Upper bounds 
We begin with the following well-known fact. 
Proposition 2. Suppose that 

(3.1) Ee zX < exp(-B$(z)), 

for some function $ : [0, oo) — > [0, oo), B > and all z > 0. Then 

(3.2) P(A > x) < e -** {x \ 
where $* = §* B is defined by 

(3.3) = sup{zx-B$(z) : z > 0}. 

Note that if <3> is an Orlicz function (a convex, continuous, non-decreasing function, such 
that $(0) =0 and $(£) — > oo as t — > oo) then $* is just a complementary function to 
Proof: This is well-known; by the usual exponentiation and Markov's inequality we have 

¥(X > x) = F(e zX > e zx ) < e - zx Ee zX < e - zx e mz) = e - {zx - B ^ z)) . 



Since the right-hand side may be minimized over z we obtain (3.2) as required. □ 



One can obtain a bound on the moment generating function of R using the fact that it 



is a limit in distribution of the iterative procedure (1.2) and verifying (3.1) for every R n . 
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In the case Q n = q (1.2) takes the form 

(3.4) R n = q + M n R n _ u 

where M n is a copy of M independent of R n -i- To argue inductively, suppose that for 
some B > 



(3.5) 



Ee 



zR n - 



1 < exp(B$(z)), z>0. 



Then by (3.4) and (3.5) applied conditionally on M n we have 



Ee 



D<?2 Ee zMn " Rn_1 < e qz Ke B ®( zMnS> 



The inductive step will be complete once we show that 

In terms of the distribution /i of M, the above inequality reads 
(3.6) 



f e B ^ zt ^{dt)<e B ^ z \ 
Jo 



Once this inequality is established, the induction is complete as one can start with arbitrary 



random variable Ro, so in particular we can ensure that (3.5) holds for R . The above 
inequality is crucial for establishing the upper bound. 

We will be interested in the tail bounds for large values of x. We assume that $ is 
non-degenerate ($(£) ^ for t ^ 0) and satisfies $(t)/t — > oo as t — > oo (i.e. $ is an iV- 
function in the language of [15J). Then $* has the same properties and it follows directly 
from the definition (3.3) that as x — > oo the supremum in (3.3) is attained at z — > oo. This 



means that it suffices that (3.1) and thus (3.6) hold only for large values of z. Thus, we 



have the following consequence of the above discussion: 



Proposition 3. Let R be given by (1.1) with Q = q. Suppose that there exist B > and 
z such that (3.6) is satisfied for the distribution of M for all z > Zq. Then 

\n¥(R > x) 



(3.7) 



lim sup 



< -1. 



4. Beta distributions 



As earlier we will denote by \x the distribution of M. Goldie-Griibel [HI Theorem 3.1] 
showed that if Q is bounded and \i and the uniform distribution on [0, 1] are equivalent at 
1 in the sense of (2.3) then the resulting perpetuity has Poissonian tails, that is 

lnP( J R>a;) 1 
a; In a; q 

Note that uniform and beta /3(a, 1) distributions are equivalent at 1. One might reasonably 
hope that considering other values of the second parameter of the beta distribution might 
lead to a different tail behavior of R but this is not the case. As we show below any M 
whose distribution is equivalent at 1 to a measure with polynomial density at 1 leads to 
the Poissonian tails of R. 



lim 

X— ¥oo 
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Theorem 4. Let the distribution of M and the beta(a,(3) distribution be equivalent at 1. 
Assume that Q = q > 0. Then 

\n¥(R>x) (3 
xhix q 

Proof: Note that all beta distributions with the same parameter and different a param- 
eters are equivalent in the sense of (2.3). Consequently, we assume for convenience that 
a = 1 so that we consider the beta distribution with the density 

f(t) = /3(i-tf-\ < t < 1 , 

which is equivalent to the distribution of M at 1. 

We show that regardless of the value of (3 > the t ails of the resulting perpetuities 
are Poissonian. To get an upper bound we verify that (3.6) holds with = e bz for a 
suitable constant b and some B > 0. Once this is done, it follows from the discussion in 
the previous section that 



lnP(i?> x) < --lnf— ) = - -x (In x -ln( Bbe)) 
b \BbeJ b 



which implies that 
(4.1) 



lnP(i?>x) 1 
lim sup < — — . 

r^oc x In x b 



Thus we are to show that for sufficiently large z > 

(4.2) e qz [ exp{Be bzt )fi{dt) < exp(Be 

Jo 

for some positive constant B and b = q//3. To that end take an e for which (2.3) holds 
with v being a beta(l, (5) distribution. Assume a t is chosen so that t > 1 — e. We split 
the integral on the left-hand side as 

e qz [ ° exp(Be bzt )fi(dt) + e qz [ exp(Be bzt )fi(dt). 

Jo J t 

The second term, through ( 2.3[ ) is bounded by 



De qz exp{Be bz )(3 [ 

J t 



1 - tf^dt = De qz exp(Be bz )(l - t 



Pick % 
(4.3) 



o 



£0(2) > 1 — e so that 



p := De qz {l-t Q y < 1. 



In order to establish (4.2) we are to show that 

-to 



qz 



It follows from (4.3) that 



exp(Be bzt )n(dt) < (1 - p)exp(Be bz ) 
t = \-e- qz l\pjDfl^ 



8 PAWEL HITCZENKO AND JACEK WESOLOWSKI 

and thus for sufficiently large z we have that to > 1 — e. Hence, the left-hand side above, 
by (2.3) again, is bounded by 

e qz exp{Be bzt °)p{0,t ) < e qz exp{Be bzt °] 



d 

1 " ^ P e 



qz 



and we want this to be less or equal than (1 — p) exp(Be ). Divide both sides by exp(£>e' 
so that the inequality to be proved reads 



qz\ 



exp (Be bzt0 - Be bz ) (l 



D 



pe 



-qz 



<1~P. 



We drop the factor 1 — j^pe qz on the left and look at the exponent. It is 

qz + Se te(i-e-'-"(p/iW _ Be bz = qz + Be bz L-^-^im™ _ 



Set b := q/p. Since p/D < 1 we have bze~ qz ^(p/ D) 1 ^ = bze~ bz (p/ D) 1 ^ < bze~ bz < 
er 1 < In 2. Since e~ u — 1 < — u/2 for < u < In 2 we see that the expression above is 
bounded by 

£?p 1/r 



qz - Bbzp 1/l3 e bz e- bz /2 = qz[l 



2/3 



and it is clear that 



e 9Z exp (Be bzt0 - Be bz ) < exp ^(1 - ^p)) > 

can be made arbitrarily small by increasing B if necessary. In particular, we can ensure 
that it is less than 1 — p for all z not too close to 0. Thus, (4.1) is proved with b = q/{3. 

To get the matching lower bound note that using again instead of M the equivalent law 
beta(l, p) with the cdf F(t) = 1 - (1 - tf we have 

v s = l-F{l-S) = 5?. 



Thus, by (2.3) 



F(R>x) > (d 



g ln(l-c) 
? Cg\ Pln(l-cq/x) 

X J 



exp ( -P j^YZ cq/xj ^ X ~ ln ( dcq ^ 



exp ( /3 ln(1 ^ C) (a;lnx)(l + o(l)) ) . 



Hence, by letting c — > + we get that 



. AnF(R>x) ^ 
hm ml > . 

x^oo a; In a; q 



□ 
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5. Generalized beta(1,/3) distributions 



In this section we consider M's whose distributions are equivalent in the sense (2.3) to 
distribution function given by 

(5.1) F{s) = F fiiV {s) = 1 - e-^- 1 ^ 1 - 5 ))", < s < 1, 0,rj>O. 

It is elementary to verify that Fp }V is indeed a distribution function which is strictly in- 
creasing on (0,1). Furthermore, Fp^ is the distribution of a beta(l,/3) random variable 
discussed in the previous section. The family Fp tT) has the following property 



F7L = F, 



1 • 



as can be easily verified by a direct calculation. Pictures of a few such distributions with 
various parameters are given in Figures [I}j2j 

For R generated with M's with distributions equivalent to the above distribution function 
the following extension of Theorem [4] holds 

Theorem 5. Let (R n ) be given by \3.4 ) where q > and M has the distribution equivalent 
to the distribution function (5.1) for some /3,T] > 0. Let R be a limit in distribution of 
(R n ). Then 

\n¥(R>x) 

lim r = . 

x-^oo x{mxy' q 
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0.6 O.S 1.0 



0.2 



2 0.4 0.6 0.8 1.0 , 0.2 0.4 0.6 0.8 



0.6 0.8 



Figure 2. (a) The distribution Fo.2,0.1, (b) its density, (c) its inverse -F510 10 , 
and (d) its density. 



Proof: For the upper bound we will show that R satisfies Proposition [3] with <&(z) 
exp^z 1 /* 7 ) for fe's in a certain range. For this <J> we have 

<!>'■(.,.) > ( I l J^Y ~ B 



1 It] 

which can be seen by using <&* B (x) > x^o — Be te » with zq = b~ v (In x) v . It follows that 
(5.2) 



\n¥(R>x) 1 
limsup — < - — . 



To verify (3.6) we will use the same argument as before; with $(z) = exp(te 1//r? ) it becomes 

e qz f exp(Be b{zt)1/v )fx(dt) < exp(Be bzl/v ), 
Jo 

where \x is the distribution of the rv M and b and B are positive constants. Splitting the 
left-hand side, with to > 1 — e as before we have 

(5e qz [ ° exp{Be b{zt)1/v )fi{dt) + e qz [ exp{Be b{zt)1/v )(i{dt). 

JO Jtn 
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By (2.3) the second term is bounded by 

De qz exp (Be 



Choose t so that p := De qz {l - F{t )) < 1. Then 

to = F^(l-pe-* z /D) 

= F^^l-pe^/D) = l-exp{-p- 1 /\-ln(pe~ qz /D)) 1 ^) 



1 — cxp 



qz 
J 



[ """ ' _ ln(/VD)^ J/, 



and for 2 sufficiently large it follows that to > 1 — £■ Now we are to prove that 
e 9 * exp (Be** 1 " 11 *"'} /i(0, £ ) < (1 - p) exp (5e 6 



By the first part of (2.3), it is enough to show that 
(5.3) e^exp (Be^ h (V^d-^) _ i)) U 



d P ^-qz-Be^ 



D 



<1-P- 



We drop the last factor on the left-hand side as it is less that 1. For to as above z l ^{l— t l J T> ) 
is close to for z sufficiently large, so that using approximations e~ x — 1 ~ — x and then 
1 — (1 — x) v ~ x/f], both valid for x close to we see that the exponent on the left-hand 
side for z sufficiently large, is 



qz 



+ Be bzl/V ( e -^a-tl /ri ) _{\„ qz _ Bbz^e bzl/V (l - t 1 ^) 



qz-—z 1/v exp ( z lh I b 



i/v 



\np/D 

qz 



1/7/ * 



~qz- — £ 1/r? exp (z 1/r > {b - {q/P) lh }) . 

For b > (q/fl) 1 ^ the second term grows faster than linearly in z, so that as long as z is not 
too close to it can be mad e arbitrarily larger than qz. Thus, (5.3) follows. Furthermore, 
letting b — > in (5.2) we obtain that 



(5.4) 



lnP(i?>x) 

limsup — < . 

x ^oo x(lnxp q 



To get a lower bound note that, using instead of the distribution of M the equivalent cdf 
Fp tV , on noting that 



1 — Fp )T ,(l — cq/x) = exp (—/J (— \n(cq/x)) v ) = exp (— /3(lnx — ln(cg))' 7 ) 
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we get for large x 

¥{R>x) > - F^ V {1 - cg/z)))MW*) 

= exp {~ W^) mnx " ln(cg))?? + ln(d)] 

= exp ^ ln( ^~ c) x(lnx)"(l - o(l))J . 
Upon letting c — ?• + it implies that 



An¥(R>x) (3 

liminf - > 

x->oo ximx) 11 q 



Combining this with (5.4) completes the proof. □ 



6. Weilbull-like tails 

In this section we explicitly construct M's that will lead to a rather different tail behavior 
of R than discussed in the previous sections. As we will see a much more rapid variability 
of M near 1 is needed to obtain a lighter tail behavior of R. More specifically, we prove 
the following theorem. 



Theorem 6. Let 1 < r < oo. Let the distribution of M be (2.3) equivalent to the distri- 
bution v with the density 

(6.1) f v (t) OC t 1 /(r-l) e - (1 - 1 r/ l . 1 -i, )P - 1J(0il)(t)> 

Then for the perpetuity R given by (1.3) with Q = q there are constants C\,C2 such that 

p lnP(i2>a:) , \n¥(R>x) 
— oo < C\ < liminf — - — < limsup — - — < c 2 < 0. 

x^oo (x/q) r x -+oc {x/qy 

Proof: For 1 < r < oo let r* be given by 

1 1 

- + - = L 

The role of r and r* are symmetric and for notational convenience we will prove the above 
inequalities for r* = r/(r — 1) rather than r. That is to say, we will show that if M is 
equivalent to a random variable whose density is proportional to 

(6.2) t V(>-*-l) e - a-trVlrl-lJjr.-l ( f ) = f -X~ ^VCr-!) ^ ? 

then the resulting perpetuity R satisfies 

AnF(R>x) , \nF(R>x) 
— oo < Ci < liminf — - — -— — < limsup — - — -— - — < c 2 < 0, 

x^oo (x/qY x ^oo {x/q) r 

for some constants Ci, c 2 . 
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Suppose we prove that for M the condition (3.5) holds for all n > 1 with = z r and 
some B > 0. Then by elementary calculation $*(a;) = r »(^i/(r-i) > so that, 

X 



(6.3) 



P(i? > s) < exp 



r *(_g r )l/(r-l) 

and this would give the claimed behavior of the logarithm of the tail probability of R. 

To establish (3.5) via inductive argument, we need to verify that (3.6) holds in the 
present situation, that is, we want to show that for z sufficiently large 



,B{ Z ty 



fx(dt) < e 



Bz r 



o 



Take e > given by (2.3) where v has density of the form (6.2) and consider 5 £ (0,e). 



Then the left hand side of the above inequality is less than 



e qz e Bz r {l-8) r _|_ e gz 



[ e B{zty ' fl{dt) < e 1z e Bz r (l-5y + De qz f 
J 1-5 Jl-t 



Consequently, we have to show that 



(6.4) 



e qz-Bz r {l-(l-SY) + De qz-Bz- 



U{t)dt < 1 



1-5 



Note that because r > 1 and < 5 < 1, the first term can be made arbitrarily small for 
z > zo sufficiently large. We thus concentrate on the second term. The following argument 
will not only complete justification of (6.4) but will also indicate how one would be led to 
a reasonable choice of f v if it were unknown. We would want to construct a density f u on 
(0, 1) for which (6.4) holds. To this end suppose for now that the density f u were of the 
form 

f v (t) = rf- l g{f). 



Upon changing variables to s = t r the second term in (6.4) becomes 

-Bz r (l-s) 



De qz-Bz' 

Setting w = 1 — s gives 
(6.5) 

We now let 



e Bz s g(s)ds = De 



qz 



{1-5)" 



g(s)ds. 



(isy 



De 



i-(isy 



-Bz r w 



g(l — w)dw. 



g(l-w) := Ke~ 1/w \ 

where 7 is to be chosen momentarily and K = K^) is set so that 

• 1 



K- 1 



^ wJ dw. 
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Then (6.5) becomes 
(6.6) 

The integrand is 
Since the function 



KDe 



-Bz 



exp I — (Bz r w + 



w Bz r w + 



rw e- l / wJ dw. 



1 



has a minimum at (7 /( y Bz r )) 1 ^ 1+1 " > whose value is 



„ 1 



(B/)^ 1 (7f+! +7 T+ 1 ) = B~t+ 1 z r ~t 



+1 



7 + 1 



7 



7/(7+1) ' 



the quantity (6.6) is no more than 



KD exp ( zg-z r 7+i£7+i 7/(7+1) I 



„ : - 7 + 1 \ 

7^ 



which upon setting 
becomes 

It is now clear that if 
(6.7) 



7 



7 + 1 



1 i.e. 7 



1 



r — 1 



KD ^{ Z ^- Bl ' \r-l)^/r 



B = A 



r , 1 



(r - 1 



.1 — 1 



where A > 1 might depend on r, then q — B x l r rj{r — 1) ' r = q{\ — A) < 0. Therefore, for 
z > Zq we obtain further 

exp {* (g - ^ 1/r (r _ 1 ) (r _ 1)/r ) } < KDe-**-» 

Thus we conclude that for z > z the left-hand side of (6.4) is bounded by 

e - 20 (B(i-(i-«5) r )4 -1 -'?) + KDe~ zm[ - A ~ l) . 

Since the value of this expression can be made smaller than 1 by choosing zq sufficiently 
large, (6.4) follows. 

Reversing the steps, we obtain the expression for the density given in (6.2) with the 
normalizing constant K r given by 

-1 



iC 1 



exp 



1 



dv. 
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0.2 0.4 0.6 0.8 



Figure 3. The density (6.2) for (a) r = 2 and (b) r = 3. 





0.2 0.4 0.6 



Figure 4. The density (6.2) for (a) r = 1.1 and (b) r = 8. 



Finally, putting the value of B given in (6.7) into (6.3) we obtain 



F(R >x)< exp 



x 



1 



which implies that 



lim sup 



lnP( J R > x) 
(x/q) r * 



To get a lower bound for ¥(R > x) we choose 5 G (0,e) as in (2.3). Then, upon passing 
to the equivalent measure with density proportional to ( |6.2 ) we have 



ps > dK r I f 1 exp , 

! 1 _ 5 v \ (1 - £r)l/(r-l) 



dt. 



Changing variables to v = (1 — t r ) ^ yields 

Ps > K 



-dv, 



(l-(l-5)'')" 1/(r - 1) V 
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for some constant K whose value is irrelevant. Since for large vq, J Vq ^~rdv is comparable 
to e~ v ° /v^ we get, up to an unimportant constant 



(1 - (1 - 5) r ) r/(r ~ 1} 



exp 



1 



(1 - (1 - 5)*-)V(r-l) 

as the lower bound for p$. Hence, up to unimportant additive terms 

lnp 5 > ln(l - (1 -S) r " 1 ' 



r -l V V (1 - (1 - 5)r)l/(r-l) (1- (l-,J)r)V(r-l) : 

as the second term above is of dominant order for 8 — > 0. For small 5 we have 

1 - (1 - S) r = 1 - exp (rln(l - 5)) ~ -rln(l - 5), 
so that upon replacing 5 by cg/a; we get that, asymptotically 

i 

! f x ^ 

lliPcg/x > 



(— r ln(l — cg/x)) 1 /^ -1 ) \cgr 



Combining this with (2.2) we get that, asymptotically, 

ln(l - c) /_ 1/(r ~ 1} \ xln(l -c) (x\ 1/(r ~ 1} 



ln(l — cq/x) \ \cqr J J cq \cqr J 



x\ ln(l — c) 

' ( cr l/r)r* ' 



It follows that 



p lnP(i2>a;) C , „ ln(l - c) 

hmmt — ; — — - — > —77 — tt, where C = < 0. 

x^oo (x/q) r ~ rV('-i)' c r 

□ 

Remarks: 

(i) The maximal value of C/r 1 '^ -1 ) is obtained by setting c = Co where Co is the unique 
solution of the equation 

1 ln(l-c) 

+ r* ■ -± ^ = 0. 



1 — c c 
The uniqueness of the solution is elementary as the function 

<*0 == ^ 

approaches — oo as c — >■ 0+ or c — >■ 1_ and 

ft (c) = -c h r • 

\1 — c c 

The expression in the parentheses, upon letting y — 1/(1 — c), y > 1, becomes 

Iny / * lny 

y — r • — — = y 1 — r 



(y-i)/y V I/" 1 
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Since — * is decreasing for y > 1, approaches 1 as y — > 1+ and as ?/ — > oo we see that 
h'(c) has exactly one sign change (from positive to negative) on (0, 1) and that this change 
occurs at c such that 



1 



+ r* 



ln(l - c ) 



0. 



1 — Co Co 

While the above equation does not have in general the closed form solution for Co as a 
function of r (or r*), the asymptotic behavior of the constant C/r 1 /^ -1 ) as r goes to or 
oo can be traced down. Since 

* Co 

' "(l-c )ln(l-Co)' 

as r — ?• oo (and thus r* — > 1 + ) we must have Co — > 0+ at the rate 1 — Co ~ 1/r*. But then 
Co~l — 1/r* = 1/r and thus 

ln(l - c ) ln(l - 1/r) 

—77 — -r — » ~ — T77 — 7v~/ — , \ 77 — 7T = rlnll — 1/r) — > —1, as r — )■ 00. 
r 1 /^-!)^ r 1 /C'--i)(l/r)'"/('"-i) v 11 ' 

Similarly, if r — )• 1 + then Co — > 1_ in such a way that 1 
ln(l — c ) — ln(r*lnr*) 



c 



- 1 l/(r* lnr* 
ln(r* lnr*) 



Then 



r 1 /^-!)^* r 1 /(»'-i)(l - l/(r* lnr*)) r * 



since, as r — > 1 + 

r l/(r-l) 



1 + 



l/(r-l) 



l/(r-l) 



e, 



and 



1 






1. 



(ii) It might appear from the argument that the form of density (6.1) was just guessed. 
While it is true that originally this was the case, there is a heuristic argument which would 
suggest the same choice. We will explain this heuristics in the next section on a different 
example, but we would like to mention that following it in the present situation would 
essentially lead to density given by (6.1). 



7. Further example 

In this section we present one more example of perpetuity that will have extremely thin 
tails. Specifically, we will show 



Proposition 7. There exist densities /m for which the perpetuity defined by (3.4) satisfies: 
, , w „ , InP(.R > x) 1 

and 

7.2 VB<q hminf 1 ~ > 1 /q) . 

x-+oc B exp(x/ B) B 



18 



PAWEL HITCZENKO AND JACEK WESOLOWSKI 



Proof: We consider the case $(z) = z In z and we will show that Proposition [3] holds for all 
B > q. It will then follow that for all such B 



(7.3) 



x 



F(R >x)< exp(- J Bexp(— - 1)), 

B 



which will imply (7.1). We will then construct a density of M which, on one hand will 
guarantee (7.3) and, on the other hand, ensure that ps is sufficiently large so that the 
argument based on Proposition [I] will give (7.2). 

To carry out the details of that plan we are to construct a density Jm for which 

• i 



D Bzt\n(zt) 



f M {t)dt < e 



Bz In 2 



This is equivalent to 



qz 



f 1 e~ B ^ zlnz t Btz f M {t)dt < 1, 
Jo 



and it is enough to construct an Jm for which 



e- B{1 - t)zXnz f M {t)dt = e qz 



-Btz In z 



f M (l-t)dt<\. 



We now set /m(1 — t) = K exp(— h(t)), where h is a non-negative function and K 1 
Jq exp(—h(t))dt. The inequality to be established becomes 



(7.4) 



-MZ 



-Btz In 



z-h(t) dt < 



~ h ^dt. 



One is guided to a reasonable choice of h by the following heuristics. Suppose h is differ- 
entiate and chosen so that 

(7.5) Btzlnz + h(t) 

is minimized at its critical point t — t z e (0, 1) which thus satisfies 

(7.6) Bz\nz + h\t z ) = 0. 



Then the left-hand side of (7.4) is no more than 

exp (qz — Bt z zlnz — h(t z )) < exp (z(q — Bt z Inz)) . 

Since we must be able to make it arbitrarily negative (by increasing B if necessary) we 
should require that t z In z is about a constant, say t z = 1/ In z for z > e. Substituting this 
into (7.6) yields 

h' (l/lnz) = — zlnz, or withs = l/lnz, h'(s) = 



,1/s 



Thus we may take 



h(t) 



1 P l/s 



-ds, 
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0.85 0.90 0.95 1.00 



Figure 5. (a) The density ju and (b) its detail closer to 1. 



and we obtain 

f M (t) = Kexp ( - 



ds 



l-t 



< t < 1, where if 1 



(Note that t z is indeed the local minimum of (7.5).) A graph of the density Jm is given in 
Figure |5j 

For the lower bound, as 



ps = K f e- h ^dt = K fe-^dt, 

J 1-5 JO 



we obtain 



lnP(i2 > x) 



ln(l — cq/x 
ln(l - c 



o 

cq/x 



x In I / e-Wdt 
cq Wo 



We need the following lemma which we justify below. 
Lemma 8. 

y In (jl ,v e- h ®dt 



(7.7) 



— > —1, as y — > oo. 



Using this lemma with y = x/(cq) and c = S/g we get, asymptotically, 
\n¥(R>x) \n(l-B/q) ( f B/x , 



which implies (7.2). 



□ 
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Proof of Lemma |$p We re- write the left-hand side of (7.7) as 

In (fv /y e- h ®dt) 

eyjy ' 

and apply l'Hospital rule. The first differentiation gives 

(-l/y 2 )e- h < 1 /»> _ e~ h ^e-y/{l - y) 

fn /y e-Wdt) (-ev/y 2 + ev/y) ~~ e-Wdt 



Differentiating again we get 

(l/y 2 )h'(l/y)e- h ^e-y/(l - y) + e~ h ^f y (e"V(l - y)) 



-ti 



1\ e~ y 



■l/y2) e -A(i/») 
o d ( e- y 



y 



yj 1-y dy \1 - y / 
Since h'(s) = —e 1 ' s /s the first term goes to — 1 as y — > oo while the second is o(l). 



□ 
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